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Theoretical and phcnomenological studies indicate that the QCD cou- 
pling as{Q^) freezes in the infrared. Hadrons may then be described by 
a perturbative expansion around "Born" states bound only by a confining 
potential. A linear potential results from the QCD equations of motion 
when Gauss' law for A° is solved with F^^i^F^" ^ as boundary condi- 
tion. The O (as°) Born states are Poincare covariant and can serve as \in) 
and {out\ states of scattering amplitudes. Their Dirac-type wave functions 
include // creation/annihilation effects giving sea-like partons at low XBj- 

1. Bound states at O {a°) 

Hadrons are highly relativistic bound states. The mass difference be- 
tween excited states is of the same order as light hadron masses, which in 
turn are much larger than the n, d, s (current) quark masses. Parton dis- 
tributions reveal the relativistic motion of quarks in the nucleon, and the 
presence of a non- vanishing sea quark distribution even at low scales [2] . 

Relativistic dynamics and color confinement are often thought to im- 
ply that the QCD coupling OsiQ'^) is large at small momentum scales Q. 
Hadrons nevertheless have features which seem difficult to reconcile with a 
strongly coupled theory. To name a few: 

• Hadron spectra reflect their valence quark (qq and qqq) degrees of 
freedom. There is no firm evidence for exotic, glueball or hybrid states. 
The sea quarks do not manifest themselves in the excitation spectrum. 

• The OZI rule [3j. E.g., the (p{W20) decays predominantly to KK, 
even though this final state is barely allowed kinematically. In a strong 

* Based on work with D. D. Dietrich and M. Jarvinen [Tl. Presented at Light Cone 
2012, Krakow, Poland. 
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coupling scenario one would expect little suppression of ss -f-T- uu, dd 
transitions at the mass scale of the strange quark, ~ 100 MeV. 

• Perturbation theory explains many features of hadron production down 
to low momentum scales [4J. 

Considerations like the above motivate studying the possibility that ag is 
of moderate size even in the confinement domain. This is not as heretical as 
it may sound. Several theoretical and phenomenological studies [5] concur 
that the strong coupling freezes at a moderate value. The quark model gives 
a semi-quantitative understanding of hadrons using the perturbative QCD 
potential added to a spin-independent linear potential. Features like the 
S — A mass splitting are then explained by single gluon exchange [B] . 

How could the confining interaction be self-consistently described theo- 
retically, and combined with perturbative QCD? One possibility is to impose 
a non-vanishing boundary condition on F^^Fa'^ in the solution of Gauss' law 
[7]. This can be illustrated in QED. Taking the diagonal matrix element 
of — V^^*^(a;) = eip^'ip{x), for a state where an electron is at xi and a 
positron at X2, gives AttA^{x;xi,X2) = e/\x — xi\ — e/\x — X2\. The stan- 
dard Coulomb potential is then ^ [6^*^(031) — eA'^(a;2)] = —al\x\ — X2\. If a 
non-vanishing field strength at spatial infinity is imposed, 

lim F-,Fli\x) = (1) 

|a;|— !>oo 

the solution of Gauss' law includes a homogeneous term, 

A^{x; xi, X2) = A'^ £ ■ X + e/47T\x — xi\ — e/47r\x — X2\ (2) 

where the unit vector i{xi, X2) may depend on the positions of the electron 
and the positron but not on x. Stationarity of the action wrt. variations 
in £ sets -lE || a;i — X2- Thus arises an instantaneous confining interaction in 
the Hamiltonian [7j, 

Ha = ^ / dxdy^mt,x)\x-y\ij^i;{t,y) (3) 

Clearly we should set A = in QED. However, a similar analysis can be 
carried out for QCD, where a linear potential is called for by data, lattice 
calculations and the quark model. The boundary condition ([T]) provides a 
dimensionful parameter which is not present in the Lagrangian, and which 
can be determined to describe (although not explain) confinement. Gauge 
covariant bound states exist for color singlet qq mesons and qqq baryons. 
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The linear potential in ([s]) is of O (eA^) and thus leading compared 
to the O (e^) perturbative potential. It was found to provide a Lorentz- 
covariant framework for bound states, giving energy eigenvalues with the 
correct dependence on the CM momentum [H [8]. This is non-trivial for 
quantization at equal time, and indicates that the implementation of the 
non-vanishing boundary condition ([T]) preserves the Poincare invariance. 

The road thus appears open to a perturbative expansion. The bound 
states formed by the linear interaction potential in Q take the place of 
the free \in) and {out\ states normally used in the scattering of pointlike 
particles. In effect, one perturbatively expands around "Born terms" which 
incorporate confinement but no perturbative interactions. 



2. Wave functions 

Relativistic dynamics necessarily involves pair creation and annihilation. 
This is manifest in the sea quark distribution of the proton, which persists 
down to low scales [2]. Consequently, relativistic bound states have an 
infinite number of Fock components. This need not exclude an analytic 
description, as demonstrated by the states of an electron in a static Coulomb 
field. The bound state energies E are given by the Dirac equation, 

[ - i7° V • 7 + eA^{x) + m7°)] (j){x) = E<j){x) (4) 

which is obtained by summing all diagrams where the electron interacts 
with the external field. As was recognized early on in the "Klein paradox" 
[9], the Dirac wave function (j){x) includes e~^e~ pair effects. Time-ordering 
of the electron interactions shows that scattering into negative energy states 
corresponds to pair creation and annihilation. 

The time- independence of ^"(a;) in Q implies that the bound state 
energies E are unchanged if retarded (instead of Feynman) electron prop- 
agators are used in all diagram^ In retarded propagation only a single 
(positive or negative energy) electron is present at any time. The Dirac 
wave function (t){x) in Q describes the electron with this boundary condi- 
tion. In analogy to cross sections [10], retarded boundary conditions give 
inclusive rather than exclusive charge densities </>^(/)(a;). 

Electron pairs contribute significantly to the Dirac charge density when- 
ever the potential is strong (and the dynamics thus is relativistic). Their 
contribution generally makes the Dirac wave function unnormalizable 

^ A time-independent external field does not transmit energy. Hence the p° component 
of the electron momentum is preserved. If > — m the negative-energy pole of the 
electron propagator at p° — —yjp^ + m? is never probed. This makes the Green 
function G{p^ ,p) independent of the ie prescription at that pole [?]• 
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This holds for any potential which is a polynomial in r or in 1/r - except 
for V{r) oc 1/r. Similarly in = 1 + 1 dimensions any potential that is a 
polynomial in x or 1/x gives unnormalizable wave functions. The absence 
of the normalization condition f d^x (j)^(f>{x) = 1 makes the Dirac energy 
spectrum continuous, quite unlike the discrete Schrodinger spectrunj^ 

Fig. 1(a) shows the Dirac wave function in Q for the QED2 potential 
eA^{x) = ^e'^\x\. Since m/e = 2.5 the dynamics is nearly non-relativistic 
at low I re I, and close agreement with the corresponding solution of the 
Schrodinger equation is indeed found for e\x\ < 5 if the Dirac solution is 
normalized to unity in this region. However, the Dirac wave function starts 
to oscillate when the potential reaches twice the electron mass, e\x\ ~ 10, in- 
dicative of contributions from e^e~ pairs. Since (l){x — )• 00) ~ exp(ie^3;^/4) 
the Dirac charge density is asymptotically constant. 



Dirac cp(x) 
Schrodinger p(x) 





p(x) Schrodinger 
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m/e = 4.0 



Fig. 1. Wave functions in £) = 1 + 1 dimensions [T. (a) Comparison of the upper 
component <f{x) of the Dirac wave function Q with the Schrodinger wave function 
p{x) for m/e = 2.5. (b) Comparison of one component of the // wave function ^ 
(for m/e = 4.0) with the Schrodinger wave function (for m/e = 2.0). 

Retarded boundary conditions may plausibly be used with the instan- 
taneous linear potentiaj^ ([s]) . Then the wave function ^{x) of an // bound 
state with CM momentum P, 

IP, *) = / .0 exp [.P . + .,)m., - |0)« 

(5) 

satisfies (for mi = m2 = m.) 
iV,-{707,$(a:)}-iP-[707,$(a;)]+m[70,$(a;)] = [E-V{x)]^{x) (6) 

where V{x) = ^eA'^\x\. The wave function ^{x) is generally singular at 
E = V{x). Requiring ^>(a;) to be regular at this point makes the energy 



This important property is bypassed in most modern textbooks. See also Ref. |12| . 
But not with perturbative Coulomb photon/gluon exchange, which (for finite fermion 
masses) transmits energy as well as 3-momentum. See footnote 1. 
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spectrum discrete rather than continuous as in the Dirac case. The // 
states in Z) = 1 + 1 were found to transform correctly under boosts fl] , and 

the energy eig envalues of Q satisfy E = ^/ P"^ + [8]. Fig. 1(b) shows 
one component of the 2x2 wave function for m/e = 4, compared 

to the Schrodinger wave function with the reduced mass m/e = 2. The 
comparison is quahtatively similar to the Dirac case in Fig. 1(a). 
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Fig. 2. (a) // ground state parton distribution in D = 1 + 1 for m/e = 0.1 
(preliminary, [T]). (b) The same distribution on a logarithmic scale. The dots are 
numerical results and the curve shows an analytic approximation valid at low x bj ■ 



3. Form factors and parton distributions 

The matrix element of the electromagnetic current j^{z) = 'il)[z)^^^^[z) 
between // bound states ^ gives the form factor 

Fj^^(z) ^ {B{P,)\f{z)\A{P,)) 

= e'^^^-^^>' j da;e*(-^''^-^")-^/2T^[^y^)^/.^o^^(^)j 

Gauge invariance, dfj,F^^{z) = 0, holds as a consequence of the bound state 
equation ^ satisfied by the wave functions ^a,^b- 

The quark distribution of target state A is obtained in the Bjorken limit 
where the photon virtuality and the mass of the final state B tend to in- 
finity. Since all states have zero width (before the perturbative corrections) 
an averaging procedure needs to be applied. The relative normalization of 
the wave functions can be determined from duality between the contri- 
butions of bound states and free quarks to the imaginary parts of current 
propagators. Our preliminary result for the quark distribution of a rela- 
tivistic ground state (m/e = 0.1) in = 1 + 1 is shown in Fig. 2. The rise 
of the distribution at low xsj is attributed to // pairs, indicating again 
the inclusive nature of the wave functions obtained with retarded boundary 
conditions. 
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